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We have studied the propagation and Zener tunnehng of hght in the bi- 
nary parabohc optical waveguide array (BPOWA), which consists of two 
evanescently coupled dissimilar optical waveguides. Due to Bragg reflections, 
BPOWA attains two minibands separated by a minigap at the zone boundary. 
Various coherent superpositions of optical oscillations and Zener tunneling oc- 
cur for different parameters on the phase diagram. In particular, Bloch-Zener 
oscillation and a different type of Bloch-dipole-Zener oscillation are obtained 
by the field-evolution analysis. The results may have potential applications 
in optical splitting and waveguiding devices and shed light on the coherent 
phenomena in optical lattices. 

(c) 2013 Optical Society of America 
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1. INTRODUCTION 

Bloch oscillation (BO) [1], dipole oscillation (DO) [2], and Zener tunneling (ZT) [3] 
are three fundamental and important transport phenomena of particles or waves in 
periodic potentials of different natures (electronic, plasmonic, photonic, acoustic, mat- 
ter wave, etc.). BO was first proposed to describe the oscillatory motion of electrons 
in periodic potential under a constant force [1]. BO cannot persist forever, but decay 
through dissipation [4], defects [5], tunnehng [6], and other damping effects. Among 
them, ZT plays an important role as an interband transition [3]. If the potential is 
not linear but parabolic, the oscillation at the bottom of the parabolic band is called 
DO [2,7]. Electronic BO is hard to observe in regular lattices because the coherence 
time of electrons is shorter than the BO period, and the electro-optic method was 
proposed to observe long dephasing time of electronic BO [8] . Photonic BO is more 
easily observed due to the longer coherence length of optical wavepackets [9] . In ad- 
dition to photonic BO [9-12], there are also various analogies of electronic BO, such 
as plasmonic BO [13-15], acoustic BO [16], matter wave BO [17], etc. 

The superposition of BO and ZT, which is referred as Bloch-Zener oscillation 
(BZO), has been studied theoretically [18-20] and experimentally [21,22]. BZO have 
potential applications in optical beam splitters and interferometers [18,22], and also 
be applied to obtain different propagation patterns in two layers of optical waveguide 
ladders [23]. Steering between BO and DO has been realized in a parabohc optical 
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waveguide array [24] . To our knowledge, there has been no study on the superposition 
of DO and ZT up to now. As we know, DO occurs in the parabohc band, and ZT 
takes place in the multiband system. To connect DO and ZT, we propose to study 
a binary parabolic optical waveguide array (BPOWA), because it has two parabolic 
minibands. 

In this work, the BPOWA has two minibands with a band gap at the edge of the 
Brillouin zone. The gap can be tuned to be small enough to make sure the occurrence 
of ZT, when the difference between the onsite propagation constants of two kinds 
of waveguides are proper. Various superpositions of optical oscillations and ZT are 
shown on the phase diagram. The superpositions of BO (approximate), DO and ZT 
can be realized, which is demonstrated by the field-evolution analysis. When the 
gap increases, the tunneling rate decreases. The band structure can be demonstrated 
visually by the spatial evolution of Bloch oscillation. The combination of BO, DO 
and ZT may have potential applications in optical splitting and waveguiding devices. 

2. MODEL AND FORMULA 

The binary parabolic optical waveguide arrays (BPOWAs) are formed by two types of 
waveguides placed alternatively with spatial period 2d, as shown in Fig. 1. We consider 
the nearest-neighbor coupling with coupling constant k. The paraxial propagation 
of light beam in the BPOWA is described by the evolution equation of the modal 
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amplitudes a„ and 6„ as 
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Qniz) + K[bn-l{z) + bn{z)] = , 

bn{z) + ii[an{z) + 0^+1(2;)] = , 
where the parabohc potential profile is V{xn) — ckx^ + a, SV is the onsite propaga- 



tion constants difference of two different kinds of waveguides, and k is the nearest- 
neighbor coupling constant. By substituting the solutions a^{z) = u^e^^""^ and 
^ni^) ~ t;™e''^'"^ into Eq. (1), we obtain the eigenvalue equation in matrix form 
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where H = - KieT''^^'^ + 1)}, {«;(e'2'='^ + 1), V{xn) + ^}}. From Eq. (2), 

the dispersion relation is derived as 



'±(x, A;) = ± 'i/2k2(i + cos2A;d) + 



(3) 



There are two bands in the dispersion relation, as shown in Fig. 1(b) with a; = 0. 
When 5y = 0, there is no gap between the two bands (as shown by the dotted fines). 
While for finite difference of onsite propagation constants, e.g., SV — 0.5, a gap opens 
at the edge of the Brillouin zone. The phase diagram of BPOWA in the ^-x domain is 
shown in Fig. 2. The two minibands /3_ and /3+ are formed with four boundary lines 
/3_(a;, 0), /3_(a;,7r/2), /3+(a;, 7r/2), and /3+(a;,0). In each band, there exists a critical 



line ^_(0, 7r/2) (or ^+(0,0)), which separates DO and BO. These two critical lines 



and an additional critical line ;5_|-(0, 7r/2) separate the two bands into ten different 
regions marked from 1 to 10, which correspond to various optical oscillations. The 
correspondences between optical oscillations and phase regions are as follows: DO 
(dipole oscillation): region 2 and 4. LBO (left Bloch oscillation): region 1, 5, 7, and 9. 
RBO (right Bloch oscillation): region 3, 6, 8, and 10. The band gap is so small that ZT 
between two minibands is possible. Five different combinations of optical oscillations 
are formed: LBO-ZT-LBO (5-9), RBO-ZT-RBO (6-10), LBO-ZT-DO-ZT-RBO (1- 
2-3), LBO-RBO (7-8), and DO (4). For each combination, the oscillation of light 
beam can transit from one kind to another through Zener tunneling. In the region 
1-2-3, the superposition of BO, DO and ZT is realized, we call such superposition as 
Bloch-dipole-Zener oscillation (BDZO). For example, a typical input light beam with 
components as shown by the red profile in region 1 can undergo the BDZO. Due to 
the involvement of DO, the oscillation spatial range of BDZO is much wider than 
the usual Bloch-Zener oscillation (BZO). This property of BDZO may have potential 
application in optical waveguiding devices. 

3. Field evolution of BZO and BDZO 

The propagation of light beam for BZO and BDZO can be demonstrated by the 
field-evolution analysis. We submit an input Gaussian beam [25] 
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with transverse wavenumber ko and intensity profile 0)p, which has a discrete 
Gaussian distribution centered at the rioth waveguide with a spatial width a. The 
exponential factor exp [—iko{n — no)] captures the phase differences between the input 
beams excited at the nth and the noth waveguides. The zero phase difference ko = 
indicates that the input beam has a plane wavefront. 

The input Gaussian beam can be expanded in terms of normal modes as 

|V'(n,0))=^A^|m), (5) 

m 

where — (m|'0(n, 0)) is the constituent component of the input Gaussian beam. 
The subsequent wave function at propagation distance z is 

\^^J{n,z))^J2^rne'^"''\m). (6) 

m 

The corresponding wave function in the reciprocal space can be obtained by taking 
the following Fourier transform 

\(l>{k,z))^J^mx,z))]. (7) 

For BZO, the components of input beam and field evolution in the spatial and re- 
ciprocal space are shown in Fig. 3 for three different cases: (a)-(c) SV = 0.0, (d)-(f) 
SV — 0.5, and (g)-(i) SV — 1.0. In each case, the same input Gaussian beam is used, 
whose components hes in the region 5, which corresponds to LBO. When 6V — 0, 
that is, the onsite propagation constants of waveguides are the same, the two bands 
are merged into one. Thus the light beam undergoes LBO between the two bound- 
aries of the merged band. The contour plots of 1-0(0;, in the spatial space and 
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in the reciprocal space are shown in Figs. 3(b) and 3(c), respectively. The 
light beam oscillates in the left side of the system in the spatial space and is almost 
constantly accelerated in the reciprocal space, which are obvious features of LBO. 
When SV — 0.5, the band splits into two minibands with a small gap, which opens 
at the edge of Brillouin zone kd — 7r/2. The gap is so small that ZT is possible. As 
shown in Fig. 3(e), the input light beam undergoes LBO at the beginning. Around 
kd = 7:/2, the beam splits into two paths due to ZT. The intensities of two LBO paths 
are almost the same, and in the reciprocal space, the slanted path of momentum also 
splits into two when ZT occurs, as shown in Fig. 3(f). When 5V — 1.0, the band 
gap becomes wider, the tunneling rate becomes smaller, that is, the intensity of LBO 
path in the upper band is stronger than that in the lower band, as shown in Fig. 3(h). 
In the reciprocal space, most of the intensity of the slanted path of momentum is 
confined in the range from — 7r/2 to 7r/2, as shown in Fig. 3(i). 

As mentioned in the description of Fig. 2, if the components of input light beam 
lies in the region 1 and the band gap is properly wide, it will undergo BDZO in the 
BPOWA. Figure 4 shows the components of input light beam and the field evolution 
in the spatial and reciprocal spaces for three different cases: (a)-(c) 5V — 0.0, (d)- 
(f) SV — 0.5, and (g)-(i) SV — 1.0. When SV — 0, the two minibands emerge into 
one, and the light beam behaves like in a unitary POWA. It undergoes DO, that is, 
oscillatory motion in both the spatial and reciprocal spaces, as shown in Figs. 4(b) 
and 4(c). When SV — 0.5, the band splits into two minibands with a small gap, as 
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shown in Fig. 4(d). The gap is so small that ZT is possible. The hght beam first 
undergoes LBO in the region 1. Around kd — 7r/2, ZT occurs, and part of light 
beam tunnels to the upper band and into the region 2, where it undergoes DO. When 
the light beam reach the boundary of the parabolic confinement at the right side, 
it can also tunnel to the lower band through ZT and undergoes RBO. In this way, 
the light beam undergoes the process LBO-ZT-DO-ZT-RBO. Figures 4(e) and 4(f) 
show the superposition of LBO, ZT, and RBO in the spatial space and reciprocal 
space, respectively. By comparing Fig. 3(e) and Fig. 4(e), we can see that the spatial 
oscillation range of BDZO is much larger than that of BZO. When 5V — 1.0, the 
band gap is too large to make the obvious occurrence of ZT. The light beam just 
undergo LBO, as shown in Figs. 4(h) and 4(i). 

It is known that the band structure can be reconstructed visually through a disper- 
sion spectroscopic approach via the spatial evolution of light beam in the photonic 
lattices [21]. Figure 5(a) shows the comparison between the rescaled spatial evolution 
of LBO and the band structure [as shown by solid lines in Fig. 1(b)]. By rescaling the 
propagation distance in unit of Zso/(27r) and the waveguide index in unit of Nh/a 
and shifting the waveguide index by —1 unit, the contour plots of spatial evolution 
of LBO match with the dispersion relation curves (solid fines) . 

The trajectories of BDZO can be roughly calculated through the Hamiltonian optics 
method. From the position-dependent dispersion relation Eq. (3), the evolution of the 
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beam can be solved by using the equations of motion 



dx d(3{x, k) 



dk 
dz 



d/3{x,k) 



(8) 
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dx 



PromEq. (8), the period of BO is estimated as — 2n/\d(3/dx\ = 'n-/\ax\. Since the 
space z in OWA is analogous to the time t in quantum system, the propagation period 
of hght beam along the z direction is the same as the oscillation period. The field- 
evolution analysis results [contour plots in Fig. 5] can be fitted with the Hamiltonian 
optics results [solid line (LBO), dashed line (DO) and dotted line (RBO) in Fig. 5]. 
We need to note that these Hamiltonian optics results are calculated by ignoring 
the multipole tunneling. In fact, due to the existence of the multipole tunneling, the 
pattern of the the field evolution is complicated. The trajectories from Hamiltonian 
optics method are used to confirm the superposition of LBO, DO and RBO. However, 
it is very hard to calculate ZT by using the Hamiltonian optics method, which can not 
show the variation of intensities. We may need other methods to obtain the tunneling 
rate of ZT in the Hamiltonian optics calculation. 

4. DISCUSSION AND CONCLUSION 

Thanks to the parabohc confinements in BPOWA, the superposition of BO, DO and 
ZT is realized. Regarding the experimental realization of such BPOWA, we can resort 
to the design of binary optical superlattices [22] and parabolic optical waveguide 
arrays [24] . The two kinds of waveguides with different onsite propagation constants 
are placed alternatively. The parabolic confinement and offsets between them are 



10 



realized by careful designing the structures, which is similar to the experimental 
design in AlGaAs waveguide arrays [12]. The precise adjustment of each individual 
waveguide can be realized in nematic liquid crystals, which has a sensitive voltage- 
controlled flexibility [26]. 

As mentioned [22], in the binary optical super lattices with linear confinement, both 
synchronous and asynchronous BZO are observed. However, in BPOWA, BO is ap- 
proximate, and the period of DO depends on the position and propagation constant 
of input beam, as a result, it is almost impossible to realize synchronous BDZO. We 
can say that BDZO in BPOWA is asynchronous, except that we can make the period 
of BO Zbo and that of DO Zdo to fulfill the relation Zbo/Zbo — p/q, where p and 
q are integers. 

In conclusion, we found that the propagation and tunneling of light beam in the 
binary parabolic optical waveguide arrays have close relation with its band structure. 
The band gap between the two minibands appears at the zone edge and can be tuned 
in width through the difference of the onsite propagation constants of two types of 
waveguides. Various combinations of optical oscillations and tunneling (BO, DO and 
ZT) are formed as demonstrated on the phase diagram. For certain band gap, both 
the Bloch-Zener oscillation (BZO) and the Bloch-dipole-Zener oscillation (BDZO) are 
realized with proper input light beams, which is demonstrated by the field-evolution 
analysis. Since the spatial oscillation range of BDZO is much larger than that of BZO, 
BDZO may have potential applications in optical splitting and waveguiding devices. 
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List of Figure Captions 

Fig. 1. (Color online) (a) Schematic diagram of a binary parabolic optical waveguide 
arrays (BPOWA). (b) Dispersion relation of unitary POWA (dotted lines) and binary 
POWA (solid lines). The parameters are N = 101, a = 2.0, k = 1.0, and SV = 0.5. 
Figl.eps 

Fig. 2. (Color onhne) Phase diagram of BPOWA. There are two bands, whose bound- 
aries are /3_(x,0), ^_(a;,7r/2), /3+(x, 7r/2), and /3+{x,0). Separated by three critical 
hues /3_(0, 7r/2), /3+(0, 7r/2), and /3+(0, 0), ten regions are formed. The correspondences 
between these regions and different optical oscillations are as follows. DO (dipole os- 
cillation): 2 and 4. LEO (left Bloch oscillation): 1, 5, 7, and 9. RBO (right Bloch 
oscillation): 3, 6, 8, and 10. The parameters are N — 101, no — 20, ko — a — 1, 
a = 2.0, K = 1.0, and 5V = 0.5. Fig2.eps 

Fig. 3. (Color online) Band structures with components of input Gaussian beams, 
contour plots of \ijj{x,z)\'^ and \(f){k,z)\'^ for three cases: (a)-(c) 5V = 0.0, (d)-(f) 
SV = 0.5, and (g)-(i) SV = 1.0. The other parameters are N = 101, Uq = 20, ko = 0, 
(7 = 1, a = 2.0, K = 1.0. Fig3.eps 

Fig. 4. (Color online) Band structures with components of input Gaussian beams, 
contour plots of \ip{x,z)\'^ and \(j){k,z)\'^ for three cases: (a)-(c) 6V = 0.0, (d)-(f) 
SV — 0.5, and (g)-(i) SV — 1.0. The other parameters are the same as those in Fig. 3 
except ko — n. Fig4.eps. 
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Fig. 5. (Color online) (a) Visual demonstration of band structure through spatial 
evolution of light beam (left Bloch oscillation) in BPOWA. The parameters are the 
same as those in Fig. 3(e). (b) Comparison of the field-evolution analysis result and 
Hamiltonian optics result of BDZO. The parameters are the same as those in Fig. 4(e). 
Fig5.eps. 
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Fig. 1. (Color online) (a) Schematic diagram of a binary parabolic optical 
waveguide arrays (BPOWA). (b) Dispersion relation of unitary POWA (dotted 
lines) and binary POWA (solid lines). The parameters are N = 101, a = 2.0, 
K = 1.0, and 6V = 0.5. Figl.eps. 
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Fig. 2. (Color online) Phase diagram of BPOWA. There are two bands, whose 
boundaries are f3_{x,0), f3_{x, n/2), (3+{x,7i/2), and 0). Separated by 

three critical lines /3-(0, tt/2), /3+(0, 7r/2), and /3+(0, 0), ten regions are formed. 
The correspondences between these regions and different optical oscillations 
are as follows. DO (dipole oscillation): 2 and 4. LBO (left Bloch oscillation): 1, 
5, 7, and 9. RBO (right Bloch oscillation): 3, 6, 8, and 10. The parameters are 
= 101, no = 20, fco = TT, a = 1, a = 2.0, k = 1.0, and 6V = 0.5. Fig2.eps. 
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5V=0.0 5V=0.5 5V=1.0 




kd kd kd 



Fig. 3. (Color online) Band structures with components of input Gaussian 
beams, contour plots of and \(j){k,z)\'^ for three cases: (a)-(c) 6V = 

0.0, (d)-(f) SV = 0.5, and (g)-(i) 6V = 1.0. The other parameters are = 101, 
no = 20, fco = 0, cr = 1, a = 2.0, k = 1.0. Fig3.eps. 
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Fig. 4. (Color online) Band structures with components of input Gaussian 
beams, contour plots of \%Ij{x,z)\'^ and for three cases: (a)-(c) 5V = 

0.0, (d)-(f) SV = 0.5, and (g)-(i) SV = 1.0. The other parameters are the same 
as those in Fig. 3 except fco = vr. Fig4.eps. 
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Fig. 5. (Color online) (a) Visual demonstration of band structure through spa- 
tial evolution of light beam (left Bloch oscillation) in BPOWA. The parameters 
are the same as those in Fig. 3(e). (b) Comparison of the field-evolution anal- 
ysis result and Hamiltonian optics result of BDZO. The parameters are the 
same as those in Fig. 4(e). Fig5.eps. 
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